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We give a brief review of the physics of gluon saturation and non-linear 
QCD evolution at small values of the Bjorken-x variable. We discuss the 
ALICE capability for low-x studies at the LHC. In particular, we con- 
centrate on the heavy quark production in the CGC framework and its 
observation with the ALICE Muon Spectrometer. 

1. Heavy quark production cross section 

The heavy quark production cross section in H1-H2 collisions can be 
illustrated by the formula: 

g h u h 2 = W dx 1 dx 2 f i (x 1 ) H 'f j (x 2 ) H2 a ij (x 1 ,x 2 , s) (1) 

id 

where <Tij(xi, x 2 , s) is the perturbative partonic cross section for two hard- 
scatterings to produce a pair of heavy quarks which depends in particular 
on the available energy. This parton level cross section is calculable as a 
power series in the strong coupling a s which depends on the renormalisation 
scale. The nucleon Parton Distribution Functions (PDF) f i /j(xi/x 2 ) Hl ^ H ' 2 
gives the probability to find the parton i/j in a nucleon of type H\/H 2 as 
a function of fraction the x\jx 2 of the nucleon's momentum carried by the 
parton. The PDFs are non-perturbative objects and thus cannot be calcu- 
lated from first principles. They are extracted from the structure functions 
of the nucleon in deep inelastic scattering (DIS) experiments. The compre- 
hension of PDFs is mandatory to make any predictions about heavy flavor 
production. 

* Presented at Excited QCD 2009, winter meeting on QCD, 8-14 February 09, Za- 
kopane, Poland 

* Magdalena.Malek@cern.ch 



(1) 



2 



M.Malek 



Excited QCD, 2009 



2. Parton structure and evolution equations 

2.1. Parton distribution in linear regime 



The DGLAP (Dokshitzer-Gribov-Lipatov-Altarelli-Parisi) [TJ evolution 
equation applicable in the Bjorken limit (Q 2 — ► oo and fixed x) takes into 
account the Q 2 dependence of the PDFs. It allows a resummation of leading 
powers of [a s ln(Q 2 )] n generated by a parton cascade in a region of the 
phase space where the transverse momenta of gluons are strongly ordered: 
Q 2 » k 2 t » k 2 n _^ t » ... » fcf t . In this approach the transverse 
momenta of the initial partons are neglected. 

At high energy (or decreasing x), the probability of gluon emission increases 
as a s ln(l/x). In this regime, the transverse momenta of the partons are of 
the same order of magnitude than the longitudinal ones and thus have to be 
taken into account in the reaction dynamics description. When the gluon 
density M is not too high, the modification of the gluon distribution with 
increasing energy can be described by the BFKL (Balitsky-Fadin-Kuraev- 
Lipatov) evolution equation [2j. This equation can be used in the Regge 
limit (s — > oo and fixed Q 2 ) of QCD. The PDFs are described in terms of kt— 
unintegrated PDFs 4>(x, kt) related to the usual gluon distribution xg(x, Q 2 ) 

by xg(x,Q 2 ) ~ J® d 2 kt4>(x, k t ). This is essentially a linear equation and 
its solution at asymptotically large energies shows an exponential growth 
of the number of gluons. This unlimited growth of the gluon density leads 
ultimately to a violation of the unitarity bound for the cross section. The 
BFKL equation must be modified so as not to violate the unitarity. 

2.2. Non-linear evolution and parton saturation 

At some small enough value of x, the gluon density is so large that the 
non-linear gluon recombination becomes important. The gluon density per 
unit area of the nucleus with atomic number A is pa ~ xg^(x, Q 2 ) /A 2 ^ and 
the gluon recombination cross section is o~ gg ^ g ~ a s /Q 2 . The border line 
between the linear and non-linear regimes is given by the saturation scale 
Q s satisfying pAO~ gg —, g = 1. The gluon recombination becomes important 
when pAO- gg -> g > 1 thus Q 2 < Q 2 with Q 2 ~ a s xgji{x, Q 2 )/A 2 / 3 . Partons 
with Q 2 > Q 2 are not affected by saturation, evolution goes towards the 
dilute system, and can be described by DGLAP or BFKL equations. On the 
other hand, partons with Q 2 < Q 2 are in the saturation regime and their 
evolution is described by the non-linear JIMWLK (Jalilian Marian-Iancu- 
McLerran-Weigert-Leonidov-Kovner) equation [3]. This equation can be 
schematized by: 



dAf(k u y) 
dy 



ciM(k t , y) - a s c 2 M 2 (k t , y) 



(2) 
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with rapidity y = log(l/x) and constant parameters c\ t 2 of the order of one. 
One remarks that the gluon recombination (or quadratic term) becomes 
important when N reaches l/a s . The gluon production stops growing, 
this is the so called gluon saturation. The saturation phenomena can be 
described by an effective theory: the Color Glass Condensate (CGC). The 
CGC is a new form of nuclear matter which controls hadronic interactions at 
asymptotically large energies. The origin of the name is the following. First 
of all, the CGC is made of small- x gluons which carry the color charge. 
Next, these small- x gluons are created by partons with larger x which act 
as a frozen color source, glass, for the emission of small- x gluons. Lastly, 
condensate because the soft degrees of freedom are as densely packed as 
they can. 

3. Low-x QCD studies at LHC with ALICE 

The LHC (Large Hadron Collider) at CERN will provide, at nominal 
operating conditions, p+p, Pb+p and Pb+Pb collisions at ^/snn = 14, 
8.8 and 5.5 TeV respectively with luminosities C ~ 10 34 , 10 29 and 5xl0 26 
cm^ 2 s~ 1 . These conditions allow to access unprecedentedly low values of 
the Bjorken-x variable. The copious production of heavy quarks (charm and 
beauty) can be used as a probe of low-x QCD phenomena. At LHC ener- 
gies heavy quarks are mainly produced through gluon-gluon fusion processes 
thus their production cross sections are significantly affected by parton dy- 
namics in the small- x regime. 

ALICE (A Large Ion Collider Experiment) [I] is the experiment dedicated 
to study the properties of matter created in Pb+Pb collisions at the LHC. 
The detector contains a wide array of sub-detectors for measuring hadrons, 
leptons and photons. Table [T] shows the x regimes for the production of cc 
and bb pairs at SPS, RHIC and LHC (for ALICE) energies. As we can see 
in Tab. [T] the increase of the collision energy involves the decrease of the 
value of the Bjorken-x. 



Machine 


SPS 


RHIC 


LHC 


LHC 


System 


Pb+Pb 


Au+Au 


Pb+Pb 


P+P 


Energy (GeV) 


17 


200 


5500 


14000 


cc 


x ~ 1CT 1 


x ~ 10~ 2 


x ~ 4 x 10" 4 


x~ 2 x 10" 4 


bb 






x ~ 2 x 10~ 3 


x ~ 6 x 10 -4 



Tabic 1. Bjorken-x values for charm and beauty production at central rapidity and 
p t = at SPS, RHIC and LHC energies. Table extracted from g]. 
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3.1. Simulations: setup 

The simulation of the heavy quark production within the CGC frame- 
work is performed with the AliRoot [5] simulation code. The theoretical 
ingredients (PDF and partonic cross section) of Eq. [T] for each approach 
are discussed in next sections. The results are compared to the reference 
simulations for the ALICE experiment (PYTHIA) that are based on the 
Mangano-Nason-Ridolfi (MNR) calculation [6]. The comparison between 
p+p and Pb+p collisions allows to draw preliminary conclusions about the 
CGC observation with the ALICE muon spectrometer |7J covering the for- 
ward pseudo-rapidity region (-4< rj <-2,5). 

Collinear factorization: PYTHIA tuned to MNR 

• p+p collisions are simulated with PYTHIA tuned to reproduce heavy 
quark production cross sections obtained with MNR calculation. The 
CTEQ |8j parametrization of the proton PDF is used. 

• Pb+p collisions are simulated from p+p ones using Glauber geometric 
model [§] to calculate the number of binary collisions. Shadowing 
effects in the gluon distribution are added by hand, using a nuclear 
modification factor. This factor, not well known at LHC energies, can 
be calculated using different parametrization methods. In this study 
we used the EKS98 parametrization [10]. Within this model, the 
transverse momentum kt of partons is given by means of a Gaussian 
distribution with a zero mean value and a width (cr = | kt |) that 
depends on the considered flavor and the colliding system. For p+p 
and Pb+p collisions, this parameter is respectively fixed to 1 and 1.16 
GeV/c for charm, and 1 and 1.6 GeV/c for beauty. 

CGC approach 

In this approach the proton PDF is given by the non-saturated CTEQ 
parametrization. The lead structure is described by the PDF taking into 
account the saturation effects. More precisely, it contains the multiple scat- 
tering corrections in the partonic cross section and the non-linear terms of 
the JIMWLK equation. The initial condition for the JIMWLK equation is 
calculated for x = 10~ 2 using the McLerran-Venugoplan model |llj . 

3.2. Results on heavy quarks level 
We calculate the charm and beauty production cross section and their 
ratios for both approaches. The ratio Rpb p = N^JdN d ]dydp t defines the 
nuclear modification factor that can be calculated from experimental data. 
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approach 


dNfJdN* 


dN^/dN* 


dN^/dN™ 


dN b A p /dNf bp 


MNR with EKS98 


0.77 


0.85 


0.03 


0.04 


CGC 


0.60 


0.80 


0.10 


0.14 



Table 2. Production ratios of charm and beauty in Pb+p collisions at 
■sj s nn = 8.8 TeV normalized to the number of binary collisions N co u and cross 
section ratios between beauty and charm for rapidity y < in MNR and CGC 
approaches. 

The results are given in Tab. |2j One remarks that the values of shadowing 
factors for beauty are very close in both models. For charm, these factors 
are very different. It shows that the gluon recombination effects are more 
important in the CGC approach that in the MNR one. In fact, the cc pairs 
are produced, on average, by fusion of gluons at smaller x than the ones 

producing the bb pairs (x = J^- exp(±y n7 j)). Since small- x gluons are 
more sensitive to recombination effects, these are more visible in the charm 
case. 

Table [2] gives also the cross section ratios for beauty over charm production. 
One sees that the MNR approach always predicts a ratio smaller than 5%, 
whereas the CGC approach gives a ratio larger than 10%. 

The evolution of the nuclear modification factor Rpbp as a function of 
transverse momentum and rapidity is presented in Fig. [TJ One remarks that 




p t (GeV/c) y 

Fig. 1. Rpbp ratios as a function of p t (left part) and rapidity (right part). 

below ~4 GeV/c, in the CGC approach, the depletion for charm is much 
more important than for beauty. It can be understood by recombination 
and multiple scattering of small-x gluons which are responsible for charm 
production. This effect is less pronounced in the MNR + EKS98 approach. 
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Shadowing effects are also visible from the rapidity dependence of Rpbp- 
Both approaches show the decrease of this ratio at high rapidities but it is 
clearly more noticeable, for both flavors, in the CGC approach. 



3.3. Results on hadron and muon level 

Next, the quark fragmentation into heavy hadrons was performed. The 
decay of produced hadrons through their semi-muonic channel was applied. 
Finally, the transport of produced muons through the muon spectrometer 
was completed [12]. It is important to highlight that a full AliRoot simu- 
lation was used to transport particles and to reconstruct their tracks; thus 
a realistic signal generation from the full detector response and a complete 
clustering and tracking procedure were applied. 

Reconstructed single muon spectra are studied in this part. The shift of 
the center of mass frame with respect to the laboratory one, by 5y = 0.47 
for Pb+p collisions, is taken into account. This kinematical shift constrains 
the rapidity window where the Rpb p can be measured; the ratio will cor- 
respond to muons with y £ [-2.97,-4] for p+p collisions and muons with 
y £ [-2.5,-3.53] for Pb+p collisions. 
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Fig. 2. Rpbp ratio for muons from charm and beauty reconstructed in the ALICE 
muon spectroometer. The left (right) plot corresponds to the results obtained with 
the CGC (MNR) model. 
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The ratio Rpbp, presented in Fig. [2] as a function of pt for muons from 
charm and beauty, exhibits a similar behavior for both flavours in the two 
models. The observation of the saturation effects via the Rpbp ratio will be 
difficult at low p t (pt < 2 GeV/c) in the ALICE muon spectrometer due 
to the different experimental cuts. 
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4. Conclusions 

The main purpose of this work was to investigate the predictions of the 
CGC approach in the case of charm and beauty production rates for the 
ALICE experiment. The results were compared to the reference simulations 
of ALICE which are based on the MNR calculation. It was observed that the 
heavy quark production is affected by small- x effects in the low pt region. It 
was shown that the shadowing effects are much more significant in the case of 
the CGC approach. It was also seen that the Cronin effect is very important 
in the CGC framework, especially for charm which is more affected because 
of its smaller mass. 

The muon level results were obtained from realistic simulations containing 
the muon spectrometer response function. Finally, we conclude that to 
observe experimentally small- a; effects with the muon spectrometer of the 
ALICE detector a very challenging analysis has to be done. 
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